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ABSTRACT 



We have performed near-infrared monitoring observations of Sgr A*, the 
Galactic center radio source associated with a supermassive black hole, with the 
near-infrared camera CIAO and the 36-element adaptive optics system on the 
Subaru telescope. We observed three flares in the Ks band (2.15/im) during 220 
minutes monitoring on 2008 May 28, and confirmed the flare emission is highly 
polarized, supporting the synchrotron radiation nature of the near-infrared emis- 
sion. Clear variations in the degree and position angle of polarization were also 
detected: an increase of the degree of polarization of about 20 %, and a swing 
of the position angle of ~ 60° — 70° in the declining phase of the flares. The 
correlation between the flux and the degree of polarization can be well explained 
by the flare emission coming from hotspot(s) orbiting Sgr A*. Comparison with 
calculations in the literature gives a constraint to the inclination angle i of the 
orbit of the hotspot around Sgr A*, as 45° < i < 90° (close to edge-on). 
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Introduction 



Sagittarius A* (Sgr A*), the compact source of radio, infrared, and X-r ay emission at the 
center of our Galaxy, is associated with a su permassive black hole fSMBH; lEckart &: Genzel 



1996 



Gillessen et a 



200 



9). 



Ghez et al.ll2000l ; ISchodel et all 120021 ) weighing in at « 4 x 1O 6 M flGhez et al. 



2008 



Sin ce th e first identification of Sgr A* in the near-infrared (NIR) 



wavelengths (IGenzel et al.ll2003l ). NIR flares have been studied extensively. Such observa- 
tions have found that Sgr A* is highly variable, and periodicity of ~ 20 minutes claimed in 
monitoring observations (e.g., Genzel et al. 2003 ; Eckart et alJbood ) would place the corre- 



sponding emission region at roughly three Schwarzschild radii for the black hole mass. If this 
picture is correct, the NIR flares from Sgr A* are powerful tools that allow us to investigate 
an accretion flow and the behavior of matter and radiation in the strong gravity regime very 
close to the SMBH. 

Polarization in the NIR bands stan ds as one of the few observables of Sgr A* not exten- 
sively investigated due to its difficulties (IMeyer et al.ll2006b| ). Successful NIR polarimetry of 
Sgr A* needs (i) the occurrence of a fairly bright flare, (ii) quick measurements to trace the 
intrinsic rapid variation of polarization, and (iii) excellent and stable atmospheric conditions 
during the entire observations. Up to now, more than a dozen NIR flares have been observed 



photometrically, bu t o nly four have been observed in linear polarization (jEckart et al.l 12006 



Meyer et al. 2006a bl; lEckart et al. 2008 ). with all the observations carried out with the 
NAOS/CONICA at the Very Large Telescope (VLT). Also, the duration of the observations 
was generally not very long (about 100 minutes). Time-resolved polarized light curves carry 
specific information about the interplay between the gravitational and magnetic fields near 
the Sgr A* event horizon, because the propagation of the polarization v ector is sensitive 
to th e presence and properties of these fields along the light trajectories (IGenzel & Karas 
2007h . 



In this letter, we present the first independent polarimetric measurements using the NIR 
camera CIAO and the adaptive optics (AO) system A036 on the 8.2 m telescope Subaru^. 
This is an almost non-stop polarimetric observation of Sgr A* over 220 minutes, and we have 
detected three main flares and shorter sub-flares superimposed on the underlying main flares. 
We will show that the profile is fairly well explained in the context of a hotspot orbiting the 
SMBH. 



1 Based on data collected at Subaru Telescope, which is operated by the National Astronomical Obser- 
vatory of Japan 
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2. Observations and Data Analysis 



On 2008 May 28, we conduc ted Ks band (2.15/i m) polarime tric observations of Sgr A* 
using CIAO and it s polarimeter (jTamura et al.ll2003f) and A036 flTakami et al.ll2004j ) on the 
Subaru telescope (live et al.l 120041 ) . CIAO provides an image of a 22? 2 x 22? 2 area of sky 
with a scale of 21.7 mas pixel -1 (Fig. [Q. With the R = 13.2 mag natural guide star USNO 
0600-28577051 located about 30" from Sgr A*, and stable atmospheric condition during the 
observations, A036 provided a stable correction with seeing values between 0?17 and 0?21 
in the K$ band. CIAO employs a rotatable half wave plate installed in front of the AO and 
a fixed wire grid polarizer in the cryostat to measure linear polarization. We made 20 sec 
exposures at four wave plate angles in the sequence of 0?0, 45?0, 22?5, and 67?5E| Including 
time for readout and rotation of the wave plate, about 200 sec is required for one sequence 
of observation. We observed a dark cloud located at a few arcmin northwest from Sgr A* 
before and after the observations of Sgr A* to obtain sky measurements. 

IRAF (Image Reduction and Analysis Facility j software package was used to perform 
dark- and flat-field corrections, sky background estimation and subtraction followed by bad 
pixel corrections. A point spread funct ion was extracted via point-source fitting with the 
program StarFinder JDiolaiti et al.lboooh fr om each image. Each im age was then deconvolved 
with the Lucy- Richardson algorithm (jLucylll974l ; lRichardsonlll972l ). and beam restored with 
a Gaussian beam of FWHM corresponding to the respective wavelengths (66 mas at 2.15 
/mi). The deconvolved images show clear flares at the position of Sgr A* (Fig. [2]). 



Flux densities of Sgr A* and other point sources were obtained via aperture photometry 
in Stokes / images [I = (I Q o + J 22?5 + I^o + J 67?5 )/2]. We used DAOFIND and APPHOT 
tasks for point sources identification and photometry. The flux density of each / image 
was calibrated wi th 10 sources in the sa me field. For extinction correction, we assumed 
Ak s = 2.8 mag (jEisenhauer et al.l 120051 ) . The flux was corrected for a background flux 
density contribution F^ g , which was determined as the mean flux measured in the same 
aperture size at six different positions, which are free from contributions of individual stars. 
The photometric error at each flux of Sgr A* was estimated by fitting a power-law to the 
rms uncertainty in the flux for non-variable faint stars of similar flux densities. We found a 



2 In this sequence of observation, rapid variation of Sgr A* flux can bias polarization measurements, but 
we confirmed that such "instrumental polarization" is only a few % except when the flux changes rapidly at 
150 minutes in the current observation (Fig. [3|). 

3 IRAF is distributed by the National Optical Astronomy Observatory which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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typical dependence of the photometric error a on flux F as a = 0.06F 78 mJy. 

The Stokes parameters /, Q, and U for Sgr A* and other point sources were determined 
by aperture polarimetry as follows. DAOFIND and APPHOT tasks were used to obtain an 
intensity for each wave plate angle (J o, J 22 o 5 , J 45 o, I 67 ° 5 ) with an aperture radius of 5 pixel 
(O'.'ll, see a circle in Fig. [2]). Based on the intensities, we calculated the Stokes parameters 
as I = (I Q o + I 22 o 5 + J 45 o + J 67 o 5 )/2, Q = I Q o — J 45 o, and U = I 22 o 5 - I 6J o 5 . With these Stokes 
parameters, the degrees P and position angles 9 of polarization for stars in the same field 
were obtained. We confirmed that P = 10.6% and 9 = 19.8° for IRS 21, and 4.8% and 24° 



for mean polarization of field stars; these agree well with the results by lOtt et al.l (119991 ). 
We obtained a small mean RMS of the degrees of polarization as 0.9% for the non-variable 
faint sources in the same field, suggesting a stable atmospheric condition and a reliable 
photometry. 

To obtain intrinsic degree and position angle of polarization for Sgr A*, contributions 
of interstellar polarization and a faint nearby star S17 should be taken into account. Here 
we assume that S17 is an intrinsically unpolarized source. An observed Stokes parameter 
normalized by intensity, (^) b , is given by 

'? j ' - (r^fr-) + (f ) • w 

1 J obs V J SgrA* + hll J V 1 J ISM 

where QsgrA* is a Stokes Q parameter for Sgr A*, isgrA* and isn are Stokes I parameters 
for Sgr A* and S17, respectively, and (^) ISM is a normalized Stokes Q parameter for the 
interstellar polarization. Hence, we can obtain an intrinsic, normalized Stokes Q parameter 
for Sgr A* as 



Q 



SgrA* 



'Sgr A* 



Q\ ( Q 



I ' ok- V I J ISM 



x ( /ssrA * + Jsi7 ) . (2) 

-^SgrA* 



(j) igM (and (j) ISM ) can be calculated accurately with bright stars in the same fields, as 
w)ism = 0-032 and (y) ISM = 0.022, in the assumption that they are located near the 
Galactic Center and are intrinsically unpolarized. We derived the last term in Eq. (j2]) as 

(-^SgrA* + -7si7)/(-fsgrA*) = (-^SgrA*+S17 ~ -^bg)/(-^SgrA*+S17 ~ (-^S17 + -^bg)), where -Fs g rA*+S17 is 

the total flux density of Sgr A* and SI 7, Fb g is a background flux density co ntribution, and 



Fsi7 = 6.7 mJy (dereddened flux; an observed magnitude is 15.3 mag, see iGillessen et al. 



20091 ). With the parameters (QsgrA*/-^SgrA*) and (£/s gr A*/-fsgrA*) an d equations 



^PsgrA* = \l ( ® SgvA ' ) + ( ^ SgrA ) ,6> Sg rA* = \ arctan 



-^SgrA* / \ -^SgrA* / - 



^SgrA* \ / ( QsgrA* 



SgrA* / / V J SgrA* 



(3) 



we can obtain the intrinsic degree PsgrA* and the position angle ^sgrA* of polarization for Sgr 
A*. The position angle 9 is measured from the north and increasing counterclockwise. 
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3. Results 

The dereddened flux of Sgr A* plus S17 show a clear variation (Fig. [3], top panel, 
black line). For reference, those for the non- variable faint sources in the same field are also 
shown in the same panel (colored lines). There are three major peaks: the broadest, and the 
weakest flare from t = to 90 minutes, the strongest flare from 145 to 190 minutes, in which 
a peak is followed by a bump, and the narrowest peak with the shortest rise/decay time 
from 195 to 210 minutes. The AO system did not work at t — 131-145 minutes; therefore, 
if there was a flare around 130 minutes, we missed that. In the first flare, substructures 
might exist around 60 minutes. The rise/decay timescale of ~ 6.5 minutes in the third flare 
is consistent with the light crossing timescales for the inner part of the accretion disk, less 
than 10 Schwarzschild radii, around a 4 x 1O 6 M black hole. 

The middle and bottom panels in Fig. [3] represent time evolutions of Q/I and U/I, 
respectively, for Sgr A* plus S17 (black line) and the non-variable faint sources (colored 
lines). The colored lines show nearly constant values. The variation of Q/I for Sgr A* is 
much larger than those of the other sources. In the equatorial coordinates, the variation 
in U/I of SgrA* is smaller than Q/I. The error bars of Q/I and U/I for Sgr A* were 
determined from the mean RMS of Q/I and U/I for the non- variable faint sources of similar 
flux densities. 

We find clear variations in the intrinsic degree of polarization PsgrA* and the position 
angle ^sgrA* of Sgr A* (Fig. Hj). The polarization was often detected with not more than 2a, 
and thus we will focus our discussion on the observations where Ps gr A* was detected with 
more than 2a (thin red lines in Fig. Hj) and more than 4a (thick red lines) levels. The first 
and second flares consist of a weakly or non- polarized main flare and a highly polarized 
sub-flare. The Ps gr A* rises slowly up to ~ 22 % and probably decays sharply in the first 
flare. In the second flare, -PsgrA* rises again up to ~ 19 %. Both rises occur after the bright 
flare phase, w hen flux of Sgr A* i s decaying. These profiles are consisten t with the previous 



observations (IMeyer et al.ll2006al ; iTrippe et al.l 120071 ; lEckart et al.ll2008l ). The variation in 



#SgrA* is complex: 9 Sgr A* swings from ~ 70° to ~ 10° in the first flare. The behaviour of 
#SgrA* between the first and second flares is not clear due to large uncertainties, but #sgrA* 
swings from ~ 70° to ~ —10° in the second flare. Due to the short duration of only ~ 13 
minutes, binning of the data points makes it impossible to trace the variation of PsgrA* and 
^sgrA* in the third flare. 
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Discussion 



So far, the time evolutions of polarization of Sgr A* in the NIR bands have been ob- 



tained just for the four flares on 2004 



( Ec kart et al.ll2006l : Mever et all20 06a 



20081 ) . As pointed out by iTrippe et al. 



u ne 12, 2005 July 29, 2006 May 31, and 2007 May 15 



b ;lTrippe et al.ll2007l ; lEckart et al.ll2008l ; IZamaninasab 



(120071 ). the observed polarization parameters, in par- 



ticular the position angle, in the first three flares show "remarkable pe rmanence", and th e 



last 2007 flare has also similar parameters in spite of some differences (lEckart et al.l 120081 ). 
In the typical flare in 2006, the degree of polarization rises up to 30% - 40%, and simultane- 



Mever et al. 


2006a: 


Tripp e et al. 


2007) 



The first flare we observed shows similar properties to those detected in the previous 
observations. It shows a rise in PsgrA* and a swing in ^sgrA* at t > 50 minutes, where the 
long flare (duration of more than 90 minutes) is slowly fading away. The observed 6*s gr A* 
during the flare pe ak is ~ 60°, whic h agrees wel l with 6 ± 20° of the flare measured by 
Eckart et all tod l. and 80 ± 10° bv iMever et aD (l2006ah . 



Such flares can be explained with a hotspot model. lEckart et al.l ( 120061 ) made fitting for 
the flares in 2004 and 2005 with a two -component (hotspot plus disk) model. For the flares 
in 2005 and 2006, Meyer et al. ( 2006al lbl) adopted a hotspot plus ring model, which leads to 
reasonable fits of polarimetric light curves. In their models, the broad overall flare is caused 
by a time varying underlying disk or ri ng, and the s h orter sub-flares are due to a hotspot on 
a relativistic orbit around the SMBH. lEckart et al.l (120081 ) proposed a temporary disk with 
a short jet model, in which quasi-periodic variation is due to a hotspot on the disk. 

The second flare from 145 to 190 minutes seems to be different from other polarized 
flares. The flare survives only for ~ 30 minutes, and is not superimposed on a broader 
underlying flare. However, similarly to the previous flares, the second flare shows sub- 
structures; the highest peak is followed by a lower peak or a bump structure. The two peaks 
are separate in time by abo ut 15 minute s , which is very similar to a quasi-periodicity of 
15.5 ± 2 minutes obtained by lMeyer et al.l (j2006bl ). so these sub-structures can be explained 
by periodic orbital motions of a single hotspot. If we consider the time evolution of degree 
of polarization, however, we can find similaritie s between our results an d the expected light 
curves for a single orbital motion of a hotspot (IBroderick &: Loebl 120061 ). 



Many authors calculated light curves from an emitting bright spot on the surface of an 



accretion di sk and com oving with the disk 



Pineault 1981; Asaoka 1989; Bao 1992 



Fukud 120031 ) . Recently, IBroderick fc Loebl (120061 ) calculated light curves in infrared wave- 
lengths, including polarization, associated with a hotspot orbiting the SMBH in the Galactic 
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center. The primary feature of the light curves with a large orbital inclination angle i (close 
to the edge-on view) is a narrow and higher peak followed by a broad lower peak/bump. 
The first peak is formed by a gravitational lensing effect which is strongest when the hotspot 
is right behind the black hole. Doppler effect and beaming due to the relativistic motion of 
the hotspot in the approaching regime make the second peak/bump. The time evolutions 
of polarized flux show a double-peak profile with a higher second peak than the first one, 
or a slow-rise and sharp-decay profile (with a lower time resolution). The variation of the 
position angle is mainly influenced by the inclination angle i. 



The similarities between the second flare and the calculations by iBroderick fc Loeb 



( 120061 ) suggest that the second flare could be explained with a single orbital motion of a 
hotspot. In the light curve of the second flare from 145 to 190 minutes, we can see the 
"first peak and second peak/bump" profile. In the time evolution of degree of polarization, 
it is clearly seen that PsgrA* increases in the decay phase of the flare. The profiles of PsgrA* 
are asymmetric, showing slow rise and sharp decay. If the observed time evolution of po- 
larization comes from such a hotspot, comparison with calculations allows us to investigate 
the inclination angle i of the hotspot orbits around Sgr A*, because the time evolution 
profile of the flux and th e degree of polarization strongly depend on the inclination angle. 
Broderick fc Loebl (120061 ) showed that when i > 67?5, light curves show a bump after the 
first peak. In addition, the time evolution of polarization has a double-peak or slow-rise and 
sharp-decay profile. These are very similar to the second flare in our observations. When 
i < 45°, by contrast, the light curve has a smooth decay phase and the polarization shows a 
symmetric or slow-decay profile. These comparisons could exclude a small incli nation angle 
(near f ace-on) of % < 45° . This is consist ent with rece nt results, i > 20° by iMeyer et al. 
(l2006ah . % > 35° bv iMever et all (bo06bh . % > 50° by iMever et al.i fj2007h . and i = 70 



with which lEckart et al.l (120081 ) found a minimum reduced-^ 2 value in the modeling of the 
observed time evolution of flux and polarization. 

Although the hotspot model is a favorite model with our results, it is premat ure to use 
such a simple model to draw strong conclusions. For example, iTrippe et al.l (120071 ) proposed 



that the swing of the position angle in NIR polarization is caused by either a magnetic field 
geometry changes due to a vanishing of the accretion disk, or materials move out of the ac- 
cretion disk, perhaps into a jet. A model of the expansion of hot self-ab sorbed synchrotron 



plasma blob was also proposed fro m multi- wavelengths observations (e.g.. lYusef-Zadeh et al. 



20061 . 120081 ; iMarrone et al.l 120081 ) . This model explains the time delay betwe e n diff erent 



waveleng t hs in flare emissions, while the hotspot model does not. IMeyer et al.l (120081 ) and 



Do et al.l (120091 ) reported non-detection of a statistically significant periodicity in NIR light 
curves. Based on the mo d el to a two-tempera t ure m agnetorotational instability driven ac- 
cretion flow by iLiu et al.l (120071 ). iHuang et al.l (120081 ) showed the spectrum and frequency- 
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dependent polarization for Sgr A* with general relativistic effects, and explain the 90° flip 
of the position angle between submillimeter and NIR observations. 



Although time-resolv ed light curves have been presented by many authors (e.g.. lGoldston et al. 



20051 ; iFalanga et al.ll2007l ). no "t ime-resolved, polar i zed" light curves have been available for 



the models other than those by iBroderick fe Loebl (120061 ). Polarization provides new infor- 



mation which is extremely useful to break degeneracy of various model parameters. Our 
observations demonstrated that it is now possible to monitor the polarimetric variation of 
Sgr A* continuously up to ~ 10 hr by combining contiguous Subaru an d VLT observations , 



as done for the NIR flux density using VLT and Keck observations ( IMeyer et al.l 120081 ) . 
Simulations including polarization evolution to test the various models for Sgr A* are now 
strongly encouraged. 

This work was supported by Grant-in- Aid for the Global COE Program " The Next Gen- 
eration of Physics, Spun from Universality and Emergence" from the Ministry of Education, 
Culture, Sports, Science and Technology (MEXT) of Japan. 
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Fig. 1. — i^s-band AO image of the central 22" of our Galaxy, obtained with the NIR 
camera CIAO on the Subaru telescope. The position of Sgr A* is indicated by the cyan 
lines. The intensity scale is logarithmic, and the integration time is 20 sec. 
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Fig. 1.— Central 1'.'3 x 1'.'3 region of the GC images in the Kg band. The time from 
the beginning of the observation is 145.1 (left) and 158.4 min (right). In each image, the 
location of Sgr A* is marked by a circle, and a faint star (S17) is also located in the circle. 
The image scale is logarithmic, and the integration time is 20 x 4 sec. The images were 
sky-subtracted, flat fielded, corrected for bad pixels and cosmic rays, and deconvolved with 
the Lucy-Richardson algorithm. 
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Fig. 3. — Top: light curves of Sgr A* plus S17 flux (thick black line), and other faint stars 
in the same field (colored thin lines). Middle and bottom: time evolution of Q/I (middle) 
and U/I (bottom) for Sgr A* plus S17 flux (thick black line), and other faint stars in the 
same field (colored thin lines). 
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Fig. 4. — Time evolutions of flux and degree of polarization (top), and flux and position 
angle (bottom) for Sgr A*. The evolution of flux is shown by black lines, and those of 
degree of polarization and position angle are shown by red lines. Contributions of S17 and 
background are subtracted in flux. The position angle is measured from the north and 
increasing counterclockwise. The data points with -PsgrA* /<7p SgrA * < 2 are shown by crosses 
with broken red lines, and those with 2 < PsgrA* l&p s rA » < 4 and -PsgrA*/o"p SgrA * > 4 are 
shown by crosses with thin and thick red lines, respectively, where PsgrA* is the intrinsic 
degree of polarization of Sgr A*, and ap , is the error of PsgrA*- 



